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| Starting with the integrated emission -measure distributions Q(T) 

f ' 

of solar active regions, we determine the distribution of the parameter 
T^ax which characterises individual plasma loops. The observed Q(T) 
were determined by combining EUV and X-ray data from two separate 
experiments on ATM/Skylab. 
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I. 


INTRODUCTION 


The structure of the ’’closed" solar corona has received much attention 
recently. Following the realisation that coronal active regions consist of 
an ensemble of plasma loops (c.f. Vaiana et. al_. , 1976), great progress has 
been made in modelling stable, quiescent loops under assumptions regarding 
the local energy balance (Rosner, Tucker and Vaiana, 1978; Craig, McClymont 
and Underwood, 1978; Emslie and Machado, 1978; Vesecky, Antiochos and Under- 
wood, 1979; Levine and Pye, 1979). Conclusions drawn by the various investi- 
gators, (c.f. Pye, et al., 1978) are that in stable loops the role of mass 
flow is negligible and that the process of non-thermal energy deposition is 
directly balanced by radiative losses; thermal conduction along a loop is 
simply the means by which heat is redistributed from one place to another. 

In papers cited above, successful comparison of predictions of the 
models with values of pressure, loop length, and maximum temperature inferred 
from observations, have 'ven confidence that energy-balance models are 
appropriate and usefu' -jcriptions of physical processes in active regions. 
These models are all very similar to one another. Antiochos (1979) and 
Habbal and Rosner (1979) have explored the conditions under which coronal 
plasma loops are stable; their work places limits upon the chromospheric 
heat flux that is required for stability (Antiochos, 1979) and upon the 
nature of the non-thermal energy supply mechanism (Habbal and Rosner, 1979). 

Given that loop stability is physically realizable and that the loops 
may be modelled on basis of local energy equilibrium, it remains to describe 
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the distribution of loops in an active region. Although instrumentation is 
available to provide resolved images of loops (c.f. Foukal, 1976) 
the nested loops of active regions are generally not resolved from one 
another. It is therefore essential that methods be developed by which the 
integrated properties of active regions may be made to yield information 
about the distribution of properties of individual loops. Levine and Pye 
(1979) have undertaken to do this; their method, based upon observational 
determination of the emission measure distribution for an active region, 
revealed some of the limitations of an emission measure analysis as it might 
be applied to this problem. 

he haye carried out a similar analysis of several coronal active regions, 
and have derived the emission measure distribution for them by combining 
EUV line fluxes with broad-band X-ray fluxes to achieve a picture of Q(T) 
across the range 4.5 < log T < 6.5. With these emission measure distribu- 
tions we are able to determine the distributions of a simple loop parameter 
- the maximum temperature T - starting with simple numerical loop models. 

The present paper describes this work, which we think holds promise for 
future development. 

II. SELECTION AND TREATMENT OF OBSERVATIONS 

Our aim, like that of Levine and Pye (1979), is to develop a method by 
which the physical properties of individual plasma loops may be determined 

in unresolved active regions by an analysis of the emission measure distri- 

> 

bution Q(T) , which is defined by 

Q(T) dT = N 2 dV. (1) 

© 
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To obtain information about Q(T) across a wide range of temperature 
we have used ATM/Skylab data, and have combined broad-band X-ray images 
from the Aero space /MSFC S-0S6 telescope with EUV raster scans obtained by 
the HCO S-055 spectroheliometer. 

a. Selection of Active Regions 

Three small, non-flaring active regions were studied. All were only a 
few days old when observed. Because they were essentially free of reported 
flaring, the derived description of the active region should refer only to 
quiescent, stable loops. Details of the observations are listed in Table I* 

b. X-Ray Observations 

We have used X-ray images photographed through two filters in the S-056 
telescope: Filter 1 (12.7 pm Al, band-pass 8-16A) and Filter 2 (6.35 p Al, 
band-pass 8-22A) . Details of the instrument are discussed by Underwood et al_. , 
(1977). The calibrations and techniques of analysis are described by Under- 
wood and McKenzie (1977). 

From the latter we directly adopted (T) , the solar spectral flux at 

the telescope focal plane per unit emission measure on the sun, for the j-th 
filter. Curves of for j = 1, 2 are shown in our Figure 1. 

2 

Copies of the flight film were microdensitometered with a 1.4 arcsec 
aperture and the photographic densities transformed via photometric calibra- 
tion to Fj , the number of deposited photons/cm“ of flight film. For a single 
pixel. 
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Figure 1. 


•*$ ^ 




Values of log gG(T) for the Three EUV Lines That Were 
Used in this Study (Given by the Left-Hand Scale) and 
Values of LogJrj for the Two X-Ray Bands Used in the 
Study (Given by the Right-Hand Scale). 
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where t = exposure time, A = telescope aperture, a = pixel area. Using the 

definition of Q(T), the full image of an active region yields the flux in 

• 1 . 

n 

F. 
t 


photons sec" 


fi d * d y 


■»/* 


(T) Q(T) dT 


(3) 


In practice there was a slight background fog; we subtracted this back- 
ground from the calculated X-ray flux to obtain corrected estimates. No 
account was taken of vignetting in the telescope. 


c . EUV Observations 

EUV rasters used by us were made in the first grating position of the 
S-055 spectroheliometer (Reeves, Huber and Timothy, 1977). We used data 
from CIII X977, OVI X1032, MgX X62S. Calibrations are given by Reeves, 


Timothy, Huber and Withbroe (1977) . 

The line fluxes in ergs sec 1 are (Withbroe, 1975): 


h ■ 4 * 


1.73 x 10 


■16 


el 


f Beff / 


G(T) Q(T)dT. 


(4) 


Atomic parameters and abundances that we adopted are given in Table II 
together with the calibration factor that converts the measured count rate 
(N counts/0.041 sec) to fluxes 1^. 

Ionization equilibrium for Be-like CIII and Li-like OVI and MgX is 
sensitive to electron density (Vernazza and Raymond, 1979). In our work we 
used their ionization equilibrium calculations for an adopted p = 1. x 10 16 K 

_3 

cm , a value based upon the results of Nicolas et a^. (1979), who obtained 
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TABLE II 


DATA FOR CALCULATION OF EUV FLUXES 


Line 

A , 
el 

f 

» 

g eff 

Calibration Factor 

C III X977 

3.8 x 10" 4 

0.14 

0.60 

1.55 

0 VI X1032 

2.1 x-10“ 4 

0.13 

0.92 

1.73 


4.6 x 10* S 4 : p 4 1.8 x 10 * 6 K cm" 3 from Si III ratios in active regions 
o 

observed with high spatial resolution. Their data may accurately reflect 

conditions at foot-points of individual loops. Vernazza and Raymond (1979) 

derived values of 10** £ N £ 2 x 10** cm 3 in active regions, suggesting 

© 

p Q -values somewhat higher than those used here; their result is however 
derived from spectra (Vernazza and Reeves, 1977) of a more energetic active 
region than the ones we have examined. 

A difficulty posed by inclusion of CIII X977 in this analysis is non- 
equilibrium ionization associated with vertical mass motions through th. 
transition region (Raymond and Dupree, 1978). Upflows can increase the line's 
intensity by as much as 30%; downflows, which are commonly seen (Raymond and 
Dupree, 1978; Nicolas et al., 1979) can decrease its intensity by a factor 
of up to 5. 

In CIII a significant population of ions is in the low-lying metastable 
3 P state and in other excited states (Dupree, 1972), so that it is necessary 
to introduce a correction factor into equation (4) to account for underpopula- 
tion of the ground state. The factor, varying from unity at low temperatures 

16 «, x 

to - 0.3 at logT= 4.5 and higher (for p Q = 10 K cm ), has been incorporated 
in our calculations. 

Values of g e ££ G(T) for each ion are plotted in Fig. 1; these include 
the excitation factor for CIII X977. 

Observed fluxes were integrated within a rectangular area on the raster 
scans that included the selected active regions and a small portion of their 
surroundings. The ratios of line brightness (active/quiet) is large enough 
(Dupree et al_. , 1973; Vernazza and Reeves, 1977) that no great error is made 


by this procedure. Indeed we found that reasonable variations of integration 
boundaries changed the measured EUV fluxes by at most five percent. Our 
solutions for Q(T) are stable against greater data fluctuations than that 
(see below). 


III. CALCULATION OF EMISSION MEASURES 

Emission measure distributions Q(T) were calculated using Withbroe's 
(1975) iterative method. Table I lists the observed fluxes and the fluxes 
calculated from the (converged) Q(T) - distributions. We considered a 
solution to be converged when the trial Q(T) had changed by less than 0.5% 
at all temperatures between iterations. While the final Q(T) do not repro- 
duce the input fluxes perfectly, the agreement between input fluxes and cal- 
culated fluxes in most cases is within 1%, except for the June 9 data, for 
which the Q(T) solution remained unconverged after 75 iterations. 

Graphs of our Q(T) are shown in Figure 2. 

We tested the stability of the method, as applied by us, by randomly 
perturbing one set of input data (Dere, 1978; Levine and Pye, 1979). With 
a random number generator we varied the input fluxes by up to 10% of their 
initial value; the Q(T) determined from the disturbed data were allowed to 
converge without limiting the number of iterations. Table III lists the means 
and standard deviations that resulted from 23 computer runs, and shoivs that 
our Q(T) are surprisingly stable against data errors. 

Our June 9 results are not secure, for two reasons: (i) a part of the 

region fell outside the raster area, and we had to estimate corrections for 
that; (ii) the EUV and X-ray observations were made 12 hours apart. This is 
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Figure 2(b). Emission Measures Q(T) for the Coronal Active Regions. 
Region 12375 at 2029 UT on 6 June 1973. 
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Figure 2(c). Emission Measures Q(T) for the Coronal Active Regions. § 

Region at 0214 UT on 8 June 1973 (solid curve), at I 

1620 UT on 8 June 1973 (dashed curve), and at 0138 UT J 

on 9 June 1973 (dot-dash curve). J 


TABLE III 


STANDARD DEVIATIONS OF Q(T) AND w m (T ) 

RANDOMLY PERTURBED DATA 0214 UT 
JUNE 8, 1973 


log T 

<Q(T)> 

0 

Q(T) 

A 
' — / 
s a 

V 

a 

w 

m 

4. 6 

4.36 

(41) cafV 1 

2.56 

(40) cnfV 1 

- 

- 

4.8 

1.38 

(42) 

8.09 

(40) 

- 

- 

5.0 

4.36 

(42) 

2.56 

(41) 

- 

- 

5.2 

1.31 

(43) 

7.72 

(41) 

- 

- 

5.4 

5.34 

(41) 

4.03 

(40) 

- 

- 

5.6 

1.09 

(41) 

1.09 

(40) 

10.1 

2.6 

5.8 

6.87 

(41) 

9.04 

(40) 

27.4 

4.1 

6.0 

1.62 

(42) 

3.54 

(41) 

14.2 

3.1 

6.2 

1.85 

(41) 

7.11 

(40) 

0.32 

0.12 

6.4 

2.12 

(40) 

3.51 

(39) 

0.0074 

0.0012 

6.6 

2.48 

(39) 

4.80 

(38) 

- 

- 

6.8 

2.75 

(38) 

9.83 

(37) 

- 

- 

7.0 

3.33 

(37) 

2.09 

(37) 

- 

- 


o 
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probably why the solution for Q(T) did not readily converge. 

Our Q(T) distributions in Figure 2 show corrugations in the temperature 
range 5.7 < log T < 6.2 that are probably artefacts of the iterative fitting 
process (Withbroe, 1975). The log T - locations of the dips correspond 

closely with regions of overlap of the log g G(T) and log /^curves of Figure 

J 

1, and likely reflect an incompatability between the EUV and X-ray calibrations 
or in the atomic data. For example, the oxygen abundance adopted by Underwood 

and McKenzie (1977) in their - calculations is a factor 3.5 lower than 

J 

the abundance adopted by us. , 

Over the range log T > 5.8 our Q(T) distributions are in reasonable agree- 
ment with the results derived by Withbroe (1975; 1977) and by Levine and Pye 
(1979), except that the maxima in our Q(T) tend to lie at lower temperatures 
then theirs. Whether this is due to the combination of EUV lines and X-ray 
bands used by us, or is an intrinsic property of such small quiescent regions 
cannot be stated. We note, however, that reported subflaring occurred only 
in region 12364, observed by us on May 30; Q(T) for this region has its maximum 
at a higher temperature than in our other distributions, in excellent accord 
with Withbroe, and Levine and Pye. 

Below log T < S.SjOur calculations show a rising Q(T) corresponding to 
plasma below the transition zone. Withbroe' s (1977) empirical emission 
measure distribution also showed an upper chromospheric rise of Q(T), but at 
lower temperatures than those obtained by us. A weakness of the present work 
that is imposed by our selection of EUV spectroheliograms for study, is that 
only one emission line has been used to probe the low temperature gas, whereas 
Withbroe used averaged values for several lines. Details of the upper 
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chromospheric Q(T) will be important for a full understanding of processes 
that take place at loop footpoints, and this is an important temperature 
range to explore. In the following discussions, we will not emphasize the 
upper chromospheric distribution of Q(T) ,oserved here. 


IV. DETERMINATION OF DISTRIBUTION OF T 

max 

One aim of an emission measure analysis is specification of the physical 
parameters that characterize the ensemble of stable loops which compose an 
unresolved active region. This approach requires that models for individual 
loops be employed that incorporate the relevant factors of distributions of 
non- thermal heating and plasma flow (if any) along a loop, the loop's geometry, 
and the heat loss rate at its footpoints. A wide range of analytic loop 
models have beon devised under various reasonable simplifying assumptions 
about these factors, and they have all yielded scaling laws of the form 
T * (p o L) B . Vesecky, Antiochos and Underwood (1979) have compared numerical 
loop models with the analytical scaling laws, and found them to be reasonably 
useful and reliable, within their range of applicability. 

Our intent is to try to develop a method by which distributions within 
active regions of one or more of the simple scaling parameters might be 
derived from an emission measure analysis. Levine and Pye undertook to do 
this with an elegant analytic model, but found they could not solve Q(T) for 
a distribution function. Like them we were not able to achieve a solution 
with the one analytic model we tried, but found that simple numerical models 
can be used. 


A ■£» 
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a. Fitting to an Analytic Model 

Emslie and Machado (1979) obtained semi-empirical loop models which 

are distinguished by two parameters, a and f, that specify the heating 

function's dependence on density (E^ - n a ) and the ratio of conductive flux 

2 4 

to radiative loss rate at the base (f = ^co^ n o ^CT)). They chose T Q = 10 K 
at the base; together with a and f this gives the spatial temperature dis- 
tribution. Finally, the emission measure q(T) is scaled by the base density 


q(T; n Q , a, f) = n Q 4> (T;a,f) (5) 

where <j> is a function given by the temperature distribution. 

We suppose that the active regions studied by us are composed of plasma 

loops whose structure may be approximated by such models, and that all loops 

are heated by the same process and have the same a. Within an active region 

the distributions of base density and ratio f are such that 

Q(T) ^ q(T; n Q , a, f) w(n o ,f) (6) 

o 

where a is fixed and w is a distribution function to be determined. We have 
no evidence for the form of w; it may be satisfactory to write 

w(n Q , f) = w n (n Q ) w f (f) (7) 


so that 


QCT) = X w f ^ (T; n o’ f) w n Cn ^j 


But, referring to equation (5), 
QCT) = 


= ]£ w £ (f) q(T; n Q , a, f)] £ w^nj) ^ 

(i J wj 


( 8 ) 


(9) 


in which the distribution of base densities merely scales the distribution 
Q(T). Thus, under assumption (7), the shape of the coronal lobe of Q(T) is 
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decided by the distribution of f. 

Although it is desirable to invert equation (9) for w.(f) this is not 
possible since q(T; n , a, f) becomes infinite at T = T . In real loops 
CVeseckey, Antiochos and Underwood, 1979) q(T) is limited by physical factors 
to finite values at T . If we suppose the semi-empirical models of Emslie 

and Machado (1979) represent loop structure just below T , we can use equation 

m 

(9) to make a rough estimate of the range of f that applies. 

5 2 

To do this we constructed q(T) using T q = 10 * K, a boundary value con- 
sistent with our observations. The highest T m are generated with a - 1 
(Emslie and Machado, 1979). For an observed Q(T) that peaks near log T = 6.0 
and in which there is still significant material at log T = 6,3, a must be 
very close to 1. 

When we adopt a = 1 we find that f = 0.92 at the peak of Q(T); since 

Q(T) is only scaled by w n Cn^) we infer from its breadth that f lies mostly 

in the range 0.89 — 0.94. This estimated range of f is in good agreement 

with Emslie and Machado (1979), even though the maximum temperatures in our 

models are lower than in theirs. Because we have used a higher boundary 
5 2 

temperature of 10 , and thus a correspondingly greater radiative loss than 

they did, our determination of the same range for f implies that a relatively 
greater conducted flux is deposited at the lower boundary. This is accomplished, 
despite the lower maximum temperatures, because less conducted flux is dissi- 
pated radiatively between the loop tops and footpoints in the models we calcu- 
lated. 

b. Fitting to Numerical Models 

In the last section we found we could not determine a distribution func- 
tion by inverting equation (9) because q(T) rises withoat limit at T . To 
investigate how this might be done we have resorted to using numerical models 

in representing q(T). 17 


We adopt a loop of constant cross-section whose emission-measure distribu- 
tion is described by its length L and maximum temperature T , 

q = q(T; L, TJ 

The observed emission measure distribution is given by 

QCT) =■ 2:2. q(T; L, T ) v (L,T 1 , (10) 

L T m m 

m 

where w is the distribution function we seek. Again we will set 


"< L -V - V«W- on 

so that 

QCT) = £ w n CV CT: L - V w L ty ] • tl2) 

We suppose that a separate analysis of magnetic fields in the active 
regions could in principle provide an insight into the character of the dis- 
tribution functions w L (L). With this information and with a model for 
q (T; L, T^) we can in principle solve equation (12) for w m ( T m )- Because we 
do not have w L (L), even in principle, we assumed it, and with adopted models 
for q(T) attempted to invert the equation. The matrix is extremely ill- 
conditioned, and we could not achieve a direct solution; instead we have 
resorted to trial -and- error estimation, a method that gives satisfactory 
results for w^ ( T m ) • In experimenting with loop models we found, however, 
that for some q(T) -models no complete solution was possible at all, even by 


this method. 


Numerical models were used which were adopted from approximate poly- 
nomial fits to the temperature and electron density distributions derived 
by Vesecky, Antiochos and Underwood (1979). With the adopted distributions 
for T(£) and N (i) we constructed q (T; L, T ) and calculated, from equation 
(12), solutions w ffl (T ) using Q(T)-data in the temperature range 5.6 < log T 
< 6.5. The calculations were restricted to this range of log T because of 
the uncertainty of Q(T) at log T < 5.5 (discussed above). Further, our sche- 
matic numerical loop models are highly simplified and are not expected to 
yield a reasonable q (T) near chromospheric temperatures. They are not neces- 
sarily in energy equilibrium everywhere, but will suffice to demonstrate 
results that can be obtained by inverting equation (12) 

From polynomial fits to the models of Vesecky, Antiochos and Underwood 
(1979) we derive approximate distributions of the form 

log (N e /N ejmax ) = -D log (V L) (13) 

log (T/T o ) = B log + C (log l) 2 (14) 

The loops must obey a scaling law of the form T^ = f (p oo> L), where 

p = p k = N k T . We used the relationship developed by Craig, McClvmont 

oo o e,max m 

and Underwood (1978) : 

T m = 3 x 10 3 (Pco L > 4/13 - ( 15 > 

Whether scaling laws of this sort are valid for stable loops must still be 

determined (Antiochos, 1979). A numerical loop model is then specified by 

adopting values for B, C, D; selection of L and T then defines N and T . 

m 6 j max o 

Finally, the loop's emission measure distribution is (A = area): 

( T - L - V - V <rl )‘ l - cis) 
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In our calculations we fixed D = 0.36. Three sets of q (T; L, T ) 
were used in attempting solution of equation (12) for w m all cal- 
culations the ranges of L and T were 5.6 < log T < 6.5 and 9 < log L < 10.3. 

m m 

Since equation (15) requires that p QQ will vary greatly across this range of 

and L, the models are not strictly compatable with the value p Q = P 0Q A 

= 10 K cm - used in deriving the observed Q(T). Some of the loop models 

from each set are shown in Figure 3. 

In the model distribution of q (T; L, T ) the curves for a given T m are 

identical in shape, and are distinguished from one another by L. For this 

reason, the bracketed quantity in equation (12) will always have the same 

T - dependence, independent of the choice of w T (L) . This, together with the 
m L 

inconsistency in p Q , is a weakness of the present numerical models that can 

be overcome by more realistic modelling. The weighting function w^fL) that 

was used by us is given in Table IV. 

The models shown in Figure 3 were used in equation (12) to solve for the 

weighting function w (T) on 0214 UT June 8; results are shown in Figure 4. 

No solution for w (T ) is possible at the lowest T when models N1 and N2 
nr m' r m 

are used, only loops with a very steep rise of q(T) at the loop tops (model 

N5) suffice to give a solution for the weighting function at the lowest T . 

The shape of the w (T ) -distribution is essentially independent of the para- 
r m m 

meters B and C which specify the numerical models studied nere, and is only 
scaled by the adopted values of w^(L). 

Finally, we obtained distributions of the weighting function w m ( T m )> 
using loop model N3 and the weighting function w^(L) from Table IV , for all 


TABLE IV 


ASSUMED WEIGHTING FUNCTION w L (L) 


log L 


9 

1 . 

9.53 

5. 

9.66 

10. 

10. 

5. 

10.33 

1 . 
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Figure 3. Emission Measure Distributions for Numerical Loop Models 
That Were Used To Solve Equation (12) for the Weighting 
Function w m (T ). Values of B and C used in equation (14) 
are: Model Nl, B =0.9, C = -0.042, Model N2, B = 0. 82, 

C = -0. 040; Model N3, B = 0. 705, C = -0. 035. The values 
of log T m are: 5.7, 5.9, 6.1, 6.3, and 6. 5. Three models 
are shown for each T m , corresponding to Log L values of 
9.0 (upper), 9.67 (middle), and 10.33 (the lowest of three 
curves) . 
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L0G T m 

Figure 4. Solutions for the Weighting Function w m (T m ) for the Different 
Loop Models. The curves for models N1 and N2 are terminated 
at values of log T below which no solution is possible. 


the data. The resulting weighting functions are shown in Figure 5. The 
stability of these solutions against errors in observed fluxes was tested 
with the Q(T) -distributions that had been calculated from randomly dis- 
turbed input data (Section III). Standard deviations of solutions for 
w m ^ ) obtained from data perturbed by up to 10% are listed in Table III, 
and indicate that the weighting functions are probably correct as to order 
of magnitude. 

These results indicate that in compact, quiescent active regions like 
those studied here, unresolved plasma loops tend to prefer maximum tempera- 
tures at or below 10^ K. Individual, resolved loops that were examined by 
Rosner, Tucker and Vaiana (1978) and by Emslie and Machado (1979) tend to 
have maximum temperatures in the vicinity of 2.5 x 10 6 K, somewhat higher 
than we determined for these unresolved regions. 

Differences in the temperature at which the peak of Q(T) is found may 
reflect the inclusion by us of fluxes derived from X-ray images, or it may 
reflect a real difference between large and small structures on the sun. 

At issue, however, is whether a useful method can be found for determining the 
distribution of loop properties in an active region from its integrated 
emission measure distribution. The present method indicates that the dis- 
tribution of maximum temperatures T looks a great deal like the coronal 
lobe of Q(T) . 


V. SUMMARY 

We have attempted to determine, from a solar active region's integrated 
emission measure distribution, something about the properties of the indivi- 
dual loops that comprise it. The present work sets some limits on such an 
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LOG T m 

Figure 5(a). Weighting Functions w m (T ) for all Observed Data. 

Region 12364 at 15 36 UT on jO May 1973 (solid curve) 
and at 1946 UT on 31 May 1973 (dashed curve). Region 
12375 at 2029 UT on 6 June 1973 (dot-dash curve). 




6 . 6.5 


LOG T m 

Figure 5(b). Weighting Functions w m (T™) for all Observed Data. 

Region 12378. See legend for Figure 2c for identifica- 
tion of curves. 


approach. 

Our principal result is that the Q(T) observed by us can be represented 
over the full range 5.6 < log T < 6.5 by the superposition of simple loop 
models, if the models incorporate a substantial rise in their individual q(T) 
near the maximum temperature T^. This in turn suggests that the unresolved 
loops may have substantial area ratios F (Vesecky, Antiochos and Underwood, 
1979), since it is this ratio that fixes the extent of the rise in q(T). 

Since the bulk of the emission measure then is contributed from the loop 

i t 

tops, the distribution of maximum temperatures has approximately the same 
shape as does the integrated Q(T). 


./ 
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